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Abstract. A new method to extract fife from high redshift intergalactic 
absorption is described, based on the distribution of Hell Lya optical 
depths in the voids in the ionization zone of quasars. A prehminary 
estimate from recent HST-STIS spectra of PKS 1935-692 at z = 3 gives 
VLbh^ = 0.0125+o;ooi4 (l-a statistical errors, for a ACDM cosmology) 
consistent with other estimates. 



1. Introduction 

In the last few years it has become possible to observe details of absorption by 
singly ionized helium. The observations combine new information about the 
history of quasars, intergalactic gas, and structure formation. 

Early observations of the Hell Lya absorption spectral region included the 
quasars Q0302-003 (z=3.285, Jackobsen et al. 1994, HS 1700+64 (z=2.72, 
Davidsen et al. 1996) and PKS 1935-692 (z=3.18 Tytler & Jackobsen 1996). 
Higher resolution (GHRS) observations of Q0302-003 Hogan, Anderson & Rugers 
1997 and HE 2347-4342 (z=2.885, Reimers et al. 1997) revealed structure in the 
absorption which could be reliably correlated with HI absorption. Heap et al. 
have followed up with STIS on Q0302-003 (1999a) and HE 2347-4342 (1999b). 
The second observation should be particularly illuminating since HE 2347-4342 
is relatively bright allowing a high resolution grating to be used. The Anderson 
et al. (1998) observations of PKS 1935-692 with STIS yield good zero level esti- 
mates important for estimating the optical depth r. Preliminary reductions of 
longer STIS integrations of PKS 1935-692 with the 0.1 A slit by Anderson et al. 
(1999) confirm the 1998 results. Taken together, these data now appear to be 
showing the cosmic ionization of helium by quasars around redshift 3. Although 
it is possible that the medium is already ionized to Helll by other sources at a 
lower level (Miralda-Escude et al. 1999). 

All of the objects show absorption with mean r > 1 at redshifts lower than 
the quasar. For the higher redshift QSO's Q0302-033 and PKS 1935-692 (shown 
in Fig. 1) there is a clear shelf of r>1.3 in a wavelength region of order 20 A 
in observed wavelength blueward of the quasar emission line redshift, dropping 
to a level consistent with zero flux or r>3 beyond that. The sharp edge led 
Anderson et al. to conclude that gas initially containing helium as mostly Hell 
was being double-ionized in a region around the quasars. The lack of a strong 
emission line for Hell Lya suggests that ionizing flux is escaping so that the 228 
A flux may be similar to a simple power-law extension of the observed 304 A 
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rest frame flux. Hogan et al. (1997) used the same reasoning to estimate the 
time required for quasars to create the double ionized helium region to be 20 
Myr for a 20 A shelf (dependent on the Hubble parameter, spectral hardness, 
cosmology, baryon density and the shelf size). 

The features present in Hell Lya spectra are reflected in the HI Lya forest 
for these quasars. Attempts to model the HeH absorption with line systems 
detected in HI suggest that very low column HI absorbers, difficult to differen- 
tiate from noise in HI spectra, provide a substantial contribution to the Hell 
absorption. Typically, in the shelf region, the ratio of HcII to HI ions is of order 
20 or more, rising to at least 100 farther away (The cross-section for Hell Lya 
absorption is 1/4 that of HI). Both PKS 1935-692 and HE 2347-4342 display 
conspicuous voids in the Hell absorption near the apparent edge of the Helll 
bubble with corresponding voids in the HI spectra. 

Rcionization and the origin of the Hell Lya forest can addressed with de- 
tailed theoretical treatments (eg. Zheng & Davidsen 1995, Zhang et al. 1998, 
Fardal et al. 1998, Abel & Haehnelt 1999, Gnedin 1999). Wadsley, Hogan & An- 
derson (1999) presented numerical models of the onset of full helium rcionization 
around a single quasar. The models integrated one-dimensional radiative trans- 
fer along lines of sight taken from cosmological hydrodynamical simulations and 
used flux levels comparable to PKS 1935-692. The results reinforced the basic 
intepretation of PKS 1935-692 by Hogan et al. (1997) as the growth of an He HI 
bubble over time in a medium that was mostly Hell. In addition void recoveries 
in the large ionized bubbles similar to that observed for PKS 1935-692 occurred 
often among a random set of simulated lines of sight. The gas temperature in 
such bubbles is strongly affected by the ionization of He II to He III, particularly 
because the first photons to reach much of the gas will be quite hard since the 
softer photons are absorbed close to the quasar until the gas is optically thin. In 
particular the underdense medium reaches temperatures of order 15000K. We 
make a simple analytical argument for this result in the next section. 

Given our fairly good guess as to the physical conditions near PKS 1935- 
692 we are in position to attempt to convert the observed optical depths into 
a density in baryons. Comparing optical depths measured for the small sample 
of voids in the shelf region of PKS 1935-692 to distributions of void widths 
and densities from simulations we can build an estimate of 0,b, the total cosmic 
density in baryons. The are still uncontrollable systematic uncertainties in this 
estimate but these differ from other techniques and can be addressed with better 
simulations and a larger sample of quasars. The main focus of this paper is the 
technique. 



2. Estimating Using Helium Quasar Observations 

In He II Lya absorption spectra only the voids are sufficiently low density to 
allow measurements of the optical depth. The highly ionized void gas is opti- 
cally thin to the ionizing radiation and cool enough to ignore collisional effects, 
resulting in a fairly simple relation between optical depth and gas density. 
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THell = 1.0 X 10^0^HeIl/(^/100 kmS-^MpC-^) 

= 106 ^{p/pf^,h^ (i±^)' (^/lOO kms-^Mpc-^)"\ (1) 



where THeii is the optical depth for absorption of He II Lya, nfjeii is the number 
density of He II ions, ^ converts from real space to wavelength (velocity) space 
and p is the gas density. The recombination rate a{T) contains the only temper- 
ature dependence. The photo-ionizations per second F is determined from the 
rest frame Hell ionizing flux, can be extrapolated from the rest frame flux 

at 304^ which is estimated (assuming a cosmological model) from the observed 
continuum flux at 304A(1 + z^), where Zq is the redshift of the quasar. 

The fraction of Hell depends on temperature roughly as T~^/^ through 
the recombination coefficient. The radiative transfer models of Wadsley et al. 
(1999) gave temperatures around 15000K for the underdense gas near quasars. 
The optical depth to a He II ionizing photon travelling 5z = 0.1 (The size of 
the PKS 1935-692 shelf region) at redshift z = 3 (SCDM, Q.^ = 0.05, Hq = 
SOkms""*^ Mpc~^) is approximately r = 62(i?^/54.4eV)^^. The optical depth 
falls to 1 for photons with energies, ~ 215 eV. Thus the Hell in the outer 
half of the shelf region will be ionized preferentially by photons in the energy 
range 100-200 eV. Injecting 100-54.4 eV per helium atom is equivalent to a 
temperature increase of 13000 K when the energy is distributed among all the 
particles. The cooling time for the underdense gas is very long, dominated by 
adiabatic expansion. Gas elsewhere will gain ~ 4000 K due to Helll reionization. 

As voids evolve they approach an attractor solution resembling empty uni- 
verses and their relative growth slows when they reach around 0.1 times the 
mean density, making that value fairly representative. We used A^-body sim- 
ulations of 3 cosmologies (Standard CDM, Open CDM and ACDM) to get 
probability distributions for the quantity we label normalized opacity O = 
{p/ / / ^ Hjj^^i^^) as a function of void width. O may be directly related 
to the void optical depths via (||) . The set of voids allows a statistical maximum 
likelihood comparison of the simulated distribution of O to the distribution of 
void widths and optical depths observed. Figure 1 shows the probability distri- 
bution of O in voids (equivalent to void optical depth) and void widths for the 
ACDM cosmology. The optical depth from the simulation was normalized so 
that mean density gas expanding with the Hubble flow gives an optical depth of 
1.0. 

The effect of gas pressure on the dynamics of gas in the voids should be 
negligible and was ignored in our simulations. Thus the gas follows the dark 
matter exactly. The simulations were performed using Hydra (Couchman 1991) 
with 64 Mpc periodic simulation volumes (~ 6000 kms~^ at z = 3) containing 
128^ particles. A preliminary convergence study indicated that the statistics of 
O are relatively insensitive to resolution however our resolution is still lower than 
that of the best current Lyman-a forest studies which also include gas. We were 
constrained by the need to include larger scales so that the large void portion of 
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Figure 1. Probability contours for void optical depths and velocity 
widths along random lines of sight in a ACDM simulation with the 5 
void optical depths and widths from PKS 1935-692 overplotted as green 
crosses. The opticals depths (in helium II Lyman-a absorption) are 
normalized so that mean density gas expanding with the Hubble flow 
at z = 3 has an optical depth of one. The probabilities are cumulative 
in velocity and are approximately indicative of the number of voids 
to be expected per unit redshift interval. The simulations (described 
in the text) predict ~ 3 voids with sizes ^ 300 kms~^ in the interval 
z = 3.09 to z = 3.15 where 5 were detected. 



our sample was reasonable. Higher resolution gives slightly larger estimates for 

The ratio of thcII to O that has the maximum likelihood gives an estimate 
for the total baryon density as follows, 



Qbh^ = 0.0178 -^/ 3.5 
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In this equation all the quantities are given as ratios with typical values estimated 
from PKS 1935-692, with p/p being the local density divided by the cosmologi- 
cal mean density, ^/^hubble local expansion rate along the line-of-sight 
compared to the mean, H{z) = Hq {Q.m{^ + + ^^cuRv(l + -2)^ + ^^a)^''^ local 




Figure 2. Data and models of PKS-1935-692 HST/STIS Data. The 
top (black) curves are ground based HI data (solid) and the model 
(dashed). The lower (red) curves are Hell data from STIS (solid) and 
the model (dashed). Non-zero optical depth in 5 voids was required to 
fit the data between z = 3.09 and z = 3.15: a region dominated by the 
Quasar UV Flux. This region should be free of associated absorbers. 



Hubble parameter value, a the spectral slope of the quasar {Fi, oc z^~") and Y 
the helium mass fraction. The value of H{z) is 400 kms~^Mpc~^ at z = 3 for 
an Hq = 50 kms~^Mpc~^ SCDM cosmology. 



3. Results for PKS 1935-692 

The exact systemic redshift of PKS 1935-692 is uncertain due to an absence of 
narrow IR emission lines. We use z = 3.19. 

We used the 5 large voids observed over the range z = 3.09 — 3.15 for 
PKS 1935-692. This range should be far enough from the quasar to avoid "asso- 
ciated" absorbers caught in the quasar outflow. The fits are shown in figure 2. 
He II and H I absorption are both plotted with the models for each overplotted 
as dashed lines. The He II spectrum was modeled with the known HI Lya ab- 
sorption lines and 5 parameter values for the optical depth in each void. The 
ratio of the HI to Hell optical depths was a sixth free parameter. The models 
were convolved with the same spectral point spread function as the observations 
and include same level of photon shot noise. The best model was determined 
using maximum likelihood and the fitting errors determined with Monte Carlo 
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realizations. The emptiest void has a reasonable probability according to the 
simulations, however since it has been suggested that it could be due to a lo- 
cal ionizing source we calculated the effect of removing it, which was a ~ 40% 
increase in the estimate for ilf,. 

The ionizing flux was inferred from the observed PKS 1935-692 continuum 
level at 304 A rest frame and extrapolated to z > 3 and 228 A. Assuming the 
quasar redshift is z = 3.19 this gives a flux of 6 x 10^'^^ er g s~^Hz~^cm~'^ (for 
standard CDM) at z = 3.125 which is in the middle of the voids. The sample 
of helium quasars is quite small and though PKS 1935-692 does not have any 
known variability all quasars are thought to be intrinsically variable due to the 
nature of the fueling and emission processes. Selection effects favour the idea 
that PKS 1935-692 is currently brighter than its average value over the last few 
thousand years (the ionization response timescale). If PKS 1935-692 has recently 
brightened then ( ^ implies that our estimate is biased upward. If there is 
additional foreground continuum absorption then our estimates are biased low. 

The results are tabulated below with 1-cr fitting errors indicated (Monte 
Carlo). Aside from the uncertainties mentioned above, there is still uncertainty 
in the temperature through the T^/^ dependence. If the gas was pre-ionized to 
Helll the underdense gas could be colder by a factor of order 2 which would 
lower the Qbh'^ estimates by 20 %. A larger sample of quasars would make the 
greatest improvement in the robustness of this measurement. 



Ho as 

kms"-"^ Mpc^-' 




A CDM 0.335 0.665 70 0.91 
Open CDM 0.37 70 0.91 
Standard CDM 1.0 50 0.67 


0125 ^^■^^^'^ 

U.Ui/O _0.0014 

n ni sfi +0.0026 

0.0186 _o.o023 

01 50 +0-0021 
U.UiOU _o.0019 



Table 1. Preliminary results for Qi^h? for 3 Cosmologies. The errors 
shown are fitting errors: systematic uncertainty is significantly greater. 
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